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GENERATING AND DECODING GRAPHICAL BAR CODES 



CROSS-REFERENCE TO RELATED APPLICATIONS 

This application relates to the following co-pending applications, each of 
which is incorporated herein by reference: U.S. Patent Application No. 09/579,070, 
filed May 25, 2000, by Doron Shaked et al., and entitled "A Method and Apparatus 
for Generating and Decoding a Visually Significant Bar Code;" U.S. Patent 
Application No. 09/429,515, filed October 28, 1999 by Renato Kresch et al., and 
entitled "System and Method for Counterfeit Protection;" U.S. Serial No. 09/728,292, 
filed December 1, 2000, by Jonathan Yen et al., and entitled "Authenticable 
Graphical Bar Codes;" U.S. Patent Application No. 09/578,843, filed May 25, 2000, 
by Doron Shaked et al., and entitled "Geometric Deformation Correction Method and 

System for Dot Pattern Images;" U.S. Serial No. , filed on even date herewith, by 

Doron Shaked et al., and entitled "Fiducial Mark Patterns for Graphical Bar Codes" 

[Attorney Docket No. 10016689-1]; and U.S. Patent Application No. , filed on 

even date herewith, by Jonathan Yen et al., and entitled "Automatically Extracting 
Graphical Bar Codes" [Attorney Docket No. 10015191-1]. 

TECHNICAL FIELD 

This invention relates to systems and methods for generating and decoding 
graphical bar codes. 

BACKGROUND 

A typical bar code symbol is a pattern of parallel bars and spaces of various 
widths that represent data elements or characters. The bars represent strings of 
binary ones and the spaces represent strings of binary zeros. A conventional "one- 
dimensional" bar code symbol contains a series of bars and spaces that vary only in 
a single dimension. One-dimensional bar code symbols have relatively small 
information storage capacities. "Two-dimensional" bar codes have been developed 
to meet the increasing need for machine-readable symbols that contain more 
information than one-dimensional bar code symbols. The information storage 
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capacity of two-dimensional bar code symbols is increased relative to one- 
dimensional bar codes by varying the bar code patterns in two dimensions. 
Common two-dimensional bar code standards include PDF417, Code 1, and 
Maxicode. One-dimensional and two-dimensional bar code symbols typically are 
5 read by optical scanning techniques (e.g., by mechanically scanned laser beams or 
by self-scanning charge-coupled devices (CCD's)) that convert a printed bar code 
symbol into electrical signals. The electrical signals are digitized and decoded to 
recover the data encoded in the printed bar code symbol. 

Bar codes may be used in a variety of applications, including low information 
10 content applications (e.g., automatic price tagging and inventory management), and 
relatively high information content applications (e.g., encoding mail addresses and 
postage for automated mail reading and mail distribution systems, and encoding 
compressed content of a printed page) . 

SUMMARY 

15 As used herein, the term "graphical bar code" broadly refers to an image that 

contains inconspicuous graphical modulations that encode embedded information. 

The invention features an inventive scheme (systems and methods) for 
generating and decoding graphical bar codes characterized by relatively high 
information capacity and pleasing visual appearances that correspond substantially 
20 to their respective base image patterns. 

In one aspect of the invention, an invertible graphical operation is applied 
between regions of a base image and information-encoding graphical templates that 
are selected from a predefined template set to produce a graphical bar code with 
regions from which graphical templates are recoverable by applying an inverse 
25 graphical operation between graphical bar code regions and corresponding base 
image regions. 

Embodiments in accordance with this aspect of the invention may include one 
or more of the following features. 

In some embodiments, the invertible graphical operation corresponds to an 
30 exclusive OR (XOR) operation. In these embodiments, XOR operations may be 
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applied between the graphical bar code regions and corresponding base image 
regions to produce the graphical templates. 

In some embodiments, each of the base image regions and the graphical 
templates has the same number of pixels. Each of the base image regions and the 
5 graphical templates preferably has the same pixel layout (e.g., a rectangular or non- 
rectangular pixel array) . 

Each graphical template preferably comprises a pattern of bright and dark 
pixels. The number of bright pixels preferably is greater than the number of dark 
pixels. Each pixel location within the predefined template set preferably has an 
10 equal probability of being a dark pixel. 

In some embodiments, the graphical templates may be ordered adaptively in 
accordance with one or more predefined rules relating to disfavored graphical 
template sequences. 

In another aspect of the invention, an invertible graphical operation is applied 
^ 15 between regions of a graphical bar code and corresponding regions of a base image 
to produce a set of measurement blocks, and information-encoding graphical 
templates corresponding to the set of measurement blocks with the highest estimated 
probability is selected from a predefined template set. 

Embodiments in accordance with this aspect of the invention may include one 
~! 20 or more of the following features. 

The invertible graphical operation preferably corresponds to an XOR 
operation. 

In some embodiments, pixel value probabilities preferably are computed for 
each of the measurement blocks. The pixel value probabilities may be computed for 
25 a given measurement block based upon a weighted average of gray value 

measurements over the given measurement block. The weighted average of gray 
values may be computed by fitting a mask (e.g., a mask having a truncated Gaussian 
profile) to the dot locations over the given measurement block. Parameters of 
probability distributions (e.g., asymmetric Laplacian distributions) that are fit to a 
30 histogram of the weighted average of gray value measurements may be estimated. 

3 



Other features and advantages of the invention will become apparent from the 
following description, including the drawings and the claims. 

DESCRIPTION OF DRAWINGS 

FIG. 1 is a block diagram of an encoding module, a decoding module, and 
5 document-handling channel through which one or more graphical bar codes may be 
transmitted. 

FIG. 2 is a flow diagram of a method of rendering a base image with a 
graphical encoding of information. 

FIG. 3A is a table mapping information to a bi-level graphical code in 
10 accordance with an XOR graphical operation that is applied between regions of a 
base image and a set of information-encoding graphical templates. 

FIG. 3B is a table mapping information to a bi-level graphical code in 
accordance with an XOR graphical operation that is applied between regions of a 
base image and an alternative set of information-encoding graphical templates. 
IS FIG. 4 is a flow diagram of a process of decoding a graphical bar code based 

upon the application of a maximum likelihood principle. 

FIG. 5 is a block diagram of a process of computing probability function 
parameters for pixel values. 

FIG. 6 is an example of a histogram computed in block 88 of FIG. 5. 
20 FIG. 7 is a flow diagram of method of estimating parameters of a pair of 

probability distribution functions fit to the histogram of FIG. 6. 

DETAILED DESCRIPTION 

In the following description, like reference numbers are used to identify like 
elements. Furthermore, the drawings are intended to illustrate major features of 
25 exemplary embodiments in a diagrammatic manner. The drawings are not intended 
to depict every feature of actual embodiments nor relative dimensions of the depicted 
elements, and are not drawn to scale. 
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Referring to FIG. 1, in one embodiment, a graphical bar code 10 may be 
generated by an encoding module 12, processed through a document-handling 
channel 14, and decoded by a decoding module 16. 

Encoding module 12 includes an encoding stage 18 and a graphic modulation 
stage 20. Encoding module 12 may be implemented as one or more program 
modules that are executable on a computer or other programmable processor. In 
operation, encoding module 12 encodes information 22 into an encoded message 24. 
For example, information 22 may be compressed in accordance with a conventional 
compression algorithm, and encoded with an error correction code. Error correction 
coding provides robustness to errors due to degradations introduced by document- 
handling channel 14. The error correction codes also may be interleaved to protect 
against burst errors. In some embodiments, encoding module 12 may be 
programmed to generate from information 22 a corroborative signed message in 
accordance with the encoding process described in U.S. Serial No. 09/728,292, filed 
December 1, 2000, by Jonathan Yen et al., and entitled "Authenticate Graphical Bar 
Codes." As explained in detail below, during the encoding stage 18, the compressed 
and error correction coded message is translated into an ordered sequence of 
information-encoding graphical templates. During the graphic modulation stage 20, 
a base image 26 is modulated in accordance with the ordered sequence of graphical 
templates to produce graphical bar code 10. In particular, an invertible graphical 
operation is applied between regions of base image 26 and the ordered sequence of 
information-encoding graphical templates to produce graphical bar code 10. Base 
image 26 may be any binary discrete graphical pattern, including a logo (e.g., a 
company logo), graphics, pictures, text, images, or any pattern that has visual 
significance. Information 22 may be embedded in the graphical design of text, 
pictures, images, borders, or the background of base image 26 to produce graphical 
bar code 10. Information 22 may be embedded in graphical bar code 10 in the form 
of a binary image (e.g., a dark and bright dot pattern), a multilevel image (e.g., a 
gray-level image), or a multilevel color image. 

In the document-handling channel 14, graphical bar code 10 may be 
converted into one or more hard copies 30 by a printing stage 28. Hard copies 30 



may be processed through a handling stage 32 before being converted into an 
electronic scanned image 36 by a scanning stage 34. Graphical bar code 10 may be 
printed by a conventional printer (e.g., a LaserJet® printer available from Hewlett- 
Packard Company of Palo Alto, California, U.S.A.) or a special-purpose label printing 
device. Hard copies 30 may be in the form of any one of a wide variety of printed 
materials, including a bank draft (or check) carrying a graphical bar code of a 
withdrawal authorization signature, a stock certificate or bond carrying a graphical 
bar code of an authenticity certification, and an envelope carrying a graphical bar 
code of postage indicia. Hard copies 30 may be scanned by a conventional desktop 
optical scanner (e.g., a ScanJet® scanner available from Hewlett-Packard Company of 
Palo Alto, California, U.S.A.) or a portable scanner (e.g., a CapShare® portable 
scanner available from Hewlett-Packard Company of Palo Alto, California, U.S.A.). 
The scanned image 36 that is acquired by the scanner and introduced into decoding 
module 16 is a degraded version of the original graphical bar code 10. These 
degradations may be generated at one or more of the document-handling channel 
stages, including the printing stage 28, the handling stage 32 (e.g., copying 
degradations, stains, folds, staples, and markings), and the scanning stage 34. 

In general, decoding module 16 includes a pre-processing stage 40, an 
alignment stage 42, a geometric correction stage 44, a probabilistic analysis stage 46, 
a graphic demodulation stage 48, and a decoding stage 50. Decoding module 16 may 
be implemented as one or more program modules that are executable on a computer 
or other programmable processor. During the pre-processing stage 40, a scanned 
graphical bar code image 52 may be located in scanned image 36, and non-bar code 
regions may be cropped and trimmed from scanned image 36. During the alignment 
stage 42, fiducial marks are detected in scanned bar code image 52. The 
configuration of the detected fiducial marks indicates the type of global deformations 
that might have been introduced into the graphical bar code during transmission of 
hard copies 30 through document handling channel 14. These global deformations 
(e.g., translational, rotational, affine and skew distortions may be corrected during 
the geometric correction stage 44, as described in U.S. Patent Application No. 
09/579,070, filed May 25, 2000, by Doron Shaked et al, and entitled "A Method and 



Apparatus for Generating and Decoding a Visually Significant Bar Code." During the 
probabilistic analysis stage 46, probability models are applied to pixel value 
measurements obtained from the aligned and geometrically corrected scanned bar 
code image 53 to produce a set of probability parameters 54. As explained in detail 
below, the probability parameters 54 are used during the graphic demodulation stage 
48 to select the most likely sequence of information-encoding graphical templates 
that corresponds to the graphical template sequence that was originally encoded into 
base image 26. The selected graphical template sequence is translated into an 
encoded message 56 that is decoded into information 22 by the decoding stage 50. 

Referring to FIGS. 2, 3A and 3B, graphical bar code 10 may be produced by 
graphically modulating base image 26 in accordance with the following graphical 
encoding scheme. Base image 26 is partitioned into a plurality of sub-matrices (or 
sub-images) (step 60). For example, if base image 26 is an M x N pixel image, it 
may be partitioned into a regular array ofOKxK pixel sub-matrices, where 0 - M 
x N/(K x K). One or more regions of base image 26 are reserved for fiducial 
markings (step 62). In one embodiment, the four corners of base image 26 may be 
reserved for fiducial markings. In other embodiments, fiducial markings may 

circumscribe base image 26, as described in U.S. Serial No. , filed on even date 

herewith, by Doron Shaked et al., and entitled "Fiducial Mark Patterns for Graphical 
Bar Codes" [Attorney Docket No. 10016689-1]. The remaining un-reserved sub- 
matrices are ordered in raster scan and converted into binary bar code matrices by 
applying an invertible graphical operation (e.g., an XOR graphical operation) 
between the un-reserved sub-matrices and a sequence of information-encoding 
graphical templates that is ordered in accordance with the information 22 to be 
encoded (step 64). In general, the invertible graphical operation should generate a 
mapping (/) from the base image sub-matrices (BI,) and the graphical templates (T.) 
to the graphical bar code regions (GBCJ that may be inverted (f ) to recover 
graphical templates (T,) from the base image sub-matrices (BIJ and the graphical bar 
code regions (GBCJ. That is, 



fiB\ 9 T) -> GBC, 



(1) 



f (BI lf GBQ -> T, 



(2) 



In one embodiment, the invertible graphical operation corresponds to an XOR 
graphical operation, which operates on a pair of input pixel values to produce an 
5 output pixel value in accordance with the following XOR function: 



First Input 
Pixel Value 


Second Input 
Pixel Value 


Output Pixel 
Value 


Dark 


Dark 


Bright 


Dark 


Bright 


Dark 


Bright 


Dark 


Dark 


Bright 


Bright 


Bright 



Table 1. XOR Graphical Operation. 



10 This XOR graphical operation satisfies the above-described invertibility feature, as 
shown in equations (3) and (4) below. 

XOR{Bl p T) -> GBQ (3) 
XOR{B\ y GBQ -> T, (4) 

15 

In accordance with this graphical encoding approach, the K pixel base image 
sub-matrices are modulated and rendered as KxK pixel graphical bar code regions. 
In other words, each base image pixel is rendered as a single dot in graphical bar 
code 10 so that the original resolution of base image 26 is preserved in graphical bar 

20 code 10. In addition, suppose we have two identical base image regions, which we 
code using two different templates. The difference between the resulting encoded 
graphical bar code regions is the same as the difference between the two 
corresponding template patterns. That is, the Hamming distance between XOR 
output patterns is identical to the distance between corresponding codes. As a result, 

25 template patterns that are sufficiently different, give rise to sufficiently different XOR 
output patterns. 
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As shown in FIG. 3A, in one graphical encoding embodiment, information 22 
may be encoded into a bi-level graphical bar code image based upon a set of 
graphical templates 66 (or halftone patterns). In this embodiment, un-reserved sub- 
matrices 68 partitioned from base image 26 are modulated by applying an XOR 
5 graphical operation between base image sub-matrices 68 and graphical templates 66. 
The graphical templates 66 are 3 x 3 pixel patterns corresponding to the thirty-six 
possible arrangements of two dark dots in a 3 x 3 pattern (only five of the thirty-six 
graphical templates are shown in FIG. 3A). In some embodiments, all thirty-six 
graphical templates may be used for encoding messages 24. In other embodiments, 
10 only thirty-two of the template patterns may be used to encode five bits. In order not 
to be biased by the selection of a predefined set of codes, the particular templates 
selected for coding may be cycled through the thirty-six templates. For example, the 
code C 36 (one of thirty-six) may be obtained from the desired code C 32 (one of thirty- 
two) in accordance with equation (5): 

15 

C 36 = mod 36 (C 32 + counter) (5) 



where the counter is incremented after every use. During decoding, C 32 is decoded 
from C 36 in accordance with equation (6) : 

20 

C 32 = mod 32 (C 36 - counter) (6) 



As shown in FIG. 3B, the template pixels need not be arranged in a regular 
rectangular array. In this embodiment, graphical templates 72 correspond to the set 

25 of thirty-six possible arrangements of either one or two dark dots in a skewed 4x2 
pixel arrangement. In general, the graphical encoding process may be extended to li- 
bit encoding mappings, where n has an integer value of 2 or greater. The 
information density of this approach increases with the number (L) of graphical 
templates in accordance with log 2 L. Thus, each of the template sets of FIGS. 3A and 

30 3B has an information density of over five bits. 
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In general, the visual quality of graphical bar code 10 improves with the 
contrast level of the graphical template patterns. Assuming the contrast between 
pure dark regions and pure bright regions is 1, the average contrast (C) between dark 
and bright graphical bar code regions 70 is given by: 



where D is the number of dark dots in each template pattern and K is the total 

10 number of dot locations in each template pattern. In the embodiments of FIGS. 3A 
and 3B, the contrast level is 5/9. In addition to constructing template sets 
characterized by relatively high contrast levels, the visual quality of graphical 
template 10 may be improved by constructing the template sets such that each pixel 
location within a given template set has an equal probability of being a dark pixel. 

15 The visual quality of graphical bar code 10 also may be improved by constructing 
templates in which the local density of dark dots is uniform. In the embodiments of 
FIGS. 3A and 3B, substantial local dark density uniformity is achieved by limiting the 
number of dark pixels in each template pattern to a relatively small fraction of the 
total number of pixels. In some embodiments, the local dark dot density uniformity 

20 may be improved by modifying the template codebook adaptively. That is, the 
graphical templates may be adaptively ordered in accordance with one or more 
predefined rules relating to disfavored graphical template sequences. For example, 
in accordance with this approach, graphical template pairs with patterns that 
increase the local dark dot density when concatenated may be avoided adaptively 

25 during the decoding process. 

Referring to FIGS. 4, 5, 6 and 7, in one embodiment, graphical bar code 10 
may be decoded as follows. 

Referring initially to FIG. 4, pixel arrays in the scanned and aligned bar code 
image are identified to individual pixels in the original bar code image 10 (step 80). 

30 Typically, the graphical bar code is scanned at a high resolution such that each 
original pixel corresponds to an array of r x r scanned pixels. Due to the visual 
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degradation that often is caused by the printing, handling and scanning processes of 
document-handling channel 14, the resulting arrays of scanned pixels are no longer 
homogenous and, instead, contain scanned pixels of different gray levels. Under 
these conditions, the borders of the scanned pixels arrays corresponding to each 
5 pixel in the original graphical bar code 10 are identified in accordance with the 

method described in U.S. Patent Application No. 09/578,843, filed May 25, 2000, by 
Doron Shaked et al, and entitled "Geometric Deformation Correction Method and 
System for Dot Pattern Images." 



10 image 52 has been identified with a single original graphical bar code pixel, the 

arrays are correlated with a mask (step 82), probability density function parameters 
are estimated (step 83), and the pixel value probability (i.e., the probability of being 
a dark (or bright) pixel) is determined for each original pixel (step 84) . In the 
following description, s is defined as an array of scanned pixels corresponding to one 

15 original pixel and the gray levels of its scanned pixels are denoted by s M where i, 



Referring to FIG. 5, a weighted average of gray values is computed for each 
pixel of each measurement block (step 86). To that end, a function of the s that 
discriminates between arrays that corresponds to dark original pixels and those 
20 corresponding to bright original pixels is applied to the scanned pixel arrays. In one 
embodiment, the linear discrimination function of equation (8) is applied to the 
scanned pixel arrays to discriminate between dark and bright pixels. 



where the coefficients c lfj are taken as grid values of the two dimensional truncated 
gaussain shaped surface of equation (9). 



As explained in detail below, after each r x r pixel array in scanned bar code 



/(*)= X 



(8) 



25 




(9) 



30 
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where x indicates the mask coordinate, c the center of the dot pattern, A the 
truncation factor, s a scalar factor, and pos(x) = x for x > 0 and 0 for x<0. The 
scalars s and A , are tuned according to the scanning resolution and the expected 
deformation strength. In one embodiment, 5 and A are set so that the mask is 0 
along a band whose width is about 1/3 dot size, and rises to half its amplitude at 
about 1/3 dot size away from the center. 

The measurements are collected, and analyzed to find a threshold indicating 
whether the original pixel corresponding to a measured pattern is dark or bright. To 
this end, a histogram of the discrimination function is computed for each pixel in the 
scanned halftoned image (step 88). The histogram measurements are a weighted 
average of gray values of aligned and geometrically corrected dot locations over a 
given image. The horizontal axis is a linear function of the weighted average of gray 
values, and the vertical axis is the histogram count axis. The right side population 
belongs to "bright" pixels and the left side population belongs to "dark" pixels. 
Based on the shape of such histograms, the distribution of f(s) may be modeled by 
bi-modal non-symmetric exponential distributions. For example, in one 
embodiment, dark (bright) pixels y = f(s) may be assumed to be distributed in 
accordance with the asymmetric Laplacian distributions of equation (10). 



The parameters of the distribution functions are estimated (step 90), Since 
the parameters cc L , a R and \x are different for the dark pixel population and the bright 
pixel population, the overall number of parametrs to be estimated in this 
probabilistic approach is six. To simplify the notations in the following description, 
dark pixels are denoted by 0 and bright pixels are denoted by 1, and the parameters 
are indexed accordingly. 

As shown in FIG. 7, in order to reduced the computational burden, the 
parameter estimation process (step 90) is divided into two steps: (1) estimating ja 0 , 
m, cc 0L and a 1R ; and (2) estimating <x 1L and ot 0R . 




(10) 
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Initially, the observations y, = f(gj are partitioned into two populations based 
upon the known expected numbers of "dark" pixels (N 0 ) and "bright" pixels (NJ 
(step 92). The "dark" population includes the smallest N 0 observations and the 
"bright" populations includes the rest (largest NJ. For each of these populations the 
5 mode (|a) is estimated by the following maximum likelihood approximated method. 
The value of \i (for the sake of simplicity we ignore the 0/1 subscript) that 
maximizes the maximum likelihood estimator for jli is computed from equation (11). 

^ = ^rg{max[£loga z " a i(l^ -^)+X lo S a * ~ a *0/ ~^)]} (H) 

ieI L i<=I R 

10 

where I L is the set of observations smaller than \i, and I R is the set of observations 
larger than ja. Since the true values of oc L and oc R are not known, they are estimated 
for each tested value of [i by equations (12) and (13) (step 94). 



ielr U^J 



1/^ = 0/^)2^-^1 (13) 

ieI R 

20 where N L = | I L | and N R = | I R | , After substituting equations (12) and (13) into 
equation (11) and simplifying the resulting expression, jT may be expressed as 
follows: 



iT - Arg{max[N L (log N L - log £ (p. - y.)) +N R (log N R - log £ (y t - jx))] } (14) 

ieI L u=I R 

25 

(Note that I L and I R - and hence N L and N R - are functions of |u.) 

The most likely value of \x is computed from equation (14) (step 98). The 
estimates for a 0L and a 1R are derived directly from equations (12) and (13) using the 
computed estimates for jw 0 and (i 1 (step 100). 
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After ja 0 , jij ,a 0L and <x 1R have been estimated, the parameters a 1L and a 0R may 
be estimated (step 102). (Hereinafter, the notations a t and a 0 are used for a 1L and 
a 0R , respectively.) All observations below |i 0 and above \i x are ignored and I M and N M 
are defined asI M ={y\\i 0 <y <\i x } and N M = |I M |. The observations {y i :ieI M } are 
5 are to be samples from a random vector with the following distribution: 



P(y i9 x t ) = n P a exp{-a y t - \i } (15) 



where x I takes the values 0 or 1 with the probabilities P 0 and P 1? respectively. The 
10 EM paradigm is used to estimate a x and cc 0 . The EM algorithm begins with initial 
values for the parameters and perform a series of updating stages that converge to a 
(local) maximum of the mean maximum likelihood objective function. 

To estimate the values of a 2 and oc 0 the conditional expectation of equation 
(16) is maximized at each stage of the algorithm. 
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E[\og Piy, ,x / )|a 0 ,a 1 ,jj= (16) 
£ P(0 ^ 0 , ai , ^0[log(P 0 ^ 



The estimation for a 1L and a 0R are the values of the estimates for a x and a 0 at the last 
20 updating stage. The EM algorithm includes an initialization step and an updating 
step. During the initialization step, for each y x the gray-level of the original pixel is 
denoted by g r 

Referring back to FIG. 4, an XOR operation is applied between the scanned 
image measurement blocks and the corresponding base image regions to recover the 
25 templates embedded in the measurement blocks (step 104). In accordance with the 
XOR modulation scheme described above, the halftone pixel corresponding to g 1 has 
probability g l to be 0 and (1- g) to be 1. The value v°(Q\y) = & is defined for each 
ieI M . During the updating step, equations (17) and (18) are computed. 
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Ey(Ob,.)U-^o) 



EoU) = i€lM v Jm . . (17) 



^(7) = '^ (18) 

2^ Ob,-) 



where p'CllyJ = l-p'tO^). Next, cc 0 (j), a,(j) and p J+1 (0\y t ) are computed from 
5 equations (19), (20) and (21), respectively. 

a 0 (j) = l/E 0 (j) (19) 
a 1 (j) = l/E,(j) (20) 

p J+l (0\y.)= />°(0 1 ^^oOQexpl-a^y)^.} (21) 

P°(0 1 y i yx 0 (j)ex V {-a 0 (j)y i } + p°(l | ^(y^xp {-a ,(./>,.} 

10 

Next, the most likely template sequence corresponding to the recovered 
template sequence is selected (step 106). Given the probability of each pixel in the 
scanned image S to originate from a dark (bright) pixel, the appropriate code words 
can be scored according to the following maximum likelihood algorithm. Let 

15 w = (w p .., w L ) be a codeword form the codebook of a given sub-image p = (p,,.., p L ) 
in I. Assume that the corresponding scanned sub-image in S has discriminative 
function values y = (y^.., y L ). The maximum likelihood score for w is given by 
s(w) = P(w|y). Since each code word has equal a priori probability, under the 
assumption that observations are stochastically independent, this score is equivalent 

20 to: 

s(w) = P(y\w) = Yl P{y, I w,.) (22) 



25 The log likelihood score is given by: 
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scr(w) = logP(y\w)=Y l logPiy i \w i )^Y i (- V ) W ' log(P(y. |0)/P(y. |1)) (23) 

i=l,L i=l,L 

where means equal up to a constant that does not depend on w. Substituting 
5 P(yjO) and P(yjl) in equation (22) for the expressions of the bi-modal non- 
symmetric exponential distribution estimated above, the score for w becomes 
scr(w) = where F(yJ is given by: 





(a ^ 
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The code word that achieves the highest score scr(w) is chosen for decoding the 
embedded message (step 106), 

Finally, the selected ordered sequence of graphical templates is translated into 
symbols in accordance with the appropriate graphical template codebook. The 
15 symbols are decoded into information 22 by the decoding stage 50 of decoding 
module 16. 

In sum, the above-described embodiments provide a novel approach for 
generating and decoding graphical bar codes that are characterized by a high 
information capacity and a visual appearance that substantially corresponds to the 
20 original base image from which the graphical bar codes are generated. 

The systems and methods described herein are not limited to any particular 
hardware or software configuration, but rather they may be implemented in any 
computing or processing environment, including in digital electronic circuitry or in 
computer hardware, firmware or software. The encoding and decoding modules 
25 may be implemented, in part, in a computer program product tangibly embodied in a 
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machine-readable storage device for execution by a computer processor. In some 
embodiments, these modules preferably are implemented in a high level procedural 
or object oriented programming language; however, the algorithms may be 
implemented in assembly or machine language, if desired. In any case, the 
programming language may be a compiled or interpreted language. The encoding 
and decoding methods described herein may be performed by a computer processor 
executing instructions organized, e.g., into program modules to carry out these 
methods by operating on input data and generating output. Suitable processors 
include, e.g., both general and special purpose microprocessors. Generally, a 
processor receives instructions and data from a read-only memory and/or a random 
access memory. Storage devices suitable for tangibly embodying computer program 
instructions include all forms of non-volatile memory, including, e.g., semiconductor 
memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic 
disks such as internal hard disks and removable disks; magneto-optical disks; and 
CD-ROM. Any of the foregoing technologies may be supplemented by or 
incorporated in specially designed ASICs (application-specific integrated circuits). 
Other embodiments are within the scope of the claims. 
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